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Abstract. The band structure of fully hydrogenated Si nanosheets and nanotubes 
arc elucidated by the use of an empirical tight-binding model. The hydrogenated Si 
sheet is a semiconductor with indirect band gap of about 2.2 eV. The symmetries of 
the wave functions allow us to explain the origin of the gap. We predict that, for 
certain chiralities, hydrogenated Si nanotubes represent a new type of semiconductor, 
one with co-existing direct and indirect gaps of exactly the same magnitude. This 
behavior is different from the Hamada rule established for non-hydrogenated carbon 
and silicon nanotubes. Comparison to an ab initio calculation is made. 
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1. Introduction 

Graphene exhibits electronic properties that conventional metals do not show, e.g., 
charge carriers with zero effective mass, anomalous quantum Hall effect, and a minimum 
conductivity when carrier concentrations tend to zero [1]. The unconventional electronic 
properties of graphene occur due to its linear energy dispersion in the vicinity of the 
Dirac point and two-dimensional structure [1, 2]. Recent theoretical studies show, 
however, that a linear energy dispersion relation is not unique to graphene [3, 4, 5, 6]. 
For instance, a flat two-dimensional silicon (Si) sheet with a honeycomb lattice (silicene) 
exhibits a linear energy dispersion too [3, 4, 5, 6]. Even buckling of the Si honeycomb 
lattice (Si tends to form sp^ bonds rather than sp^ bonds) does not remove the linearity 
in the energy dispersion relation nor does it open an energy gap [3]. More generally, 
it has been predicted that a group-IV nanosheet with a flat or buckled honeycomb 
lattice also shows a linear dispersion too [4, 5, 6, 7]. It is the underlying symmetries 
of the honeycomb lattice that are responsible for the linear energy dispersion [3]. One- 
dimensional carbon nanostructures also exhibit unconventional electronic properties, 
e.g., carbon nanotubes show different electronic behavior depending on their chirality. 
Similarly to graphene, the unconventional electronic properties are not unique to carbon 
nanotubes [3, 8]. Si nanotubes also show electronic properties that are chirality 
dependent [3, 8]. 

The theoretical studies mentioned above motivated the experimental realization of 
two-dimensional sheets other than graphene. Silicon hexagonal sheets were chemically 
exfoliated from calcium disilicide (CaSi2) [9] and Si nanoribbons have been recently 
synthesized [10, 11, 12]. There is a report of the observation of a linear electronic 
dispersion for the Si nanoribbons using angle- resolved photoemission spectroscopy [13]. 

Similarly to graphene, there is interest in modifying a silicene sheet in order to 
tailor its properties. For example, studies have shown that Si nanoribbons can display 
semiconducting and magnetic behavior [4, 14]. Hydrogenation of silicene can also open a 
gap (e.g., hydrogenated graphene (graphane) is an insulator of direct band gap [15, 16]) 
and the resulting compound should be related to polysilane [17]. There are already a 
few ab initio calculations of related structures. Indeed, density functional theory (DFT) 
calculations within the local density approximation (LDA) found that a Si sheet without 
hydrogen is a semiconductor of gap zero [7]. Upon hydrogenation, the Si sheet becomes 
an semiconductor of indirect band gap [7, 17]. As for hydrogenated Si nanotubes, a 
previous density functional tight-binding theory found that hydrogenated Si nanotubes 
are semiconductors with a gap that shows little dependence on chirality [18]. 

It is the purpose of the present work to study the band structure of hydrogenated 
Si nanosheets and nanotubes in more detail than had been done previously with the 
ab initio method. In particular, (i) previous work had focused on structures, (m) the 
calculations of the hydrogenated Si sheet were all done using LDA, which suffers from 
the band-gap problem, [in) a discussion of the wave functions was missing from the 
earlier work, and (iv) a complete study of the band structures of Si-H nanotubes has 
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not yet been carried out. The most appropriate model is the empirical tight binding 
one. 

In this paper, we will first consider a hydrogenated Si sheet (also referred as 
silicane [7]). For the hydrogenated Si sheet, we find that (i) it is an semiconductor of 
indirect band gap of about 2.2 eV; {ii) the lowest point of the conduction bands (CB) 
occurs at the M point and the highest point of the valence bands (VB) occurs at the F 
point; {in) the indirect gap is closely related to the band structure of silicene and it is 
mainly due to band-filling effect; (iv) the degeneracy at the Dirac point is not lifted by 
hydrogenating the Si sheet. We then consider the electronic structure of hydrogenated 
single- walled Si nanotubes (Si-H NTs). We find that (i) Si-H NTs of certain chiralities 
are semiconductors with co-existing direct and indirect gaps of identical magnitudes; 
(ii) the magnitude of the gaps is equal to that of the hydrogenated Si sheet. 



2. Tight-binding model of hydrogenated Si sheets 

We now describe our tight-binding model. We consider a 10 x 10 Hamiltonian with the 
following Bloch states for Si {A site). 



and similarly for the B sites and for H [C and D sites). Here a = s,Px,Py,Pz are Lowdin 
orbitals; Rj is a lattice vector of the hexagonal lattice; r is the position vector of the 
Si and H atoms in the unit cell; and is the number of unit cells. For Si atoms, we 
consider first nearest neighbors located at 

a ^ a ^ a^^„ 

xH — y — Azz, — pX + Oy —Azz, 



2V3 2-" ' V3 

- ^^x- -y - A2;z, Ox + Oy + dsmz, 

and second nearest neighbors located at 

Ox + ay +0z. Ox —ay +0z, 

ay/S „ a „ ^ ^ ay/S ^ a ^ 

^x--y + Oz, -^x+-y + Oz, 

oa/S ^ a ^ ^ a-\/3 ^ a ^ 

^x+-y + Oz, -^x--y + Oz. 

Here, Az ~ 0.19 a is the buckling height; dsm — 0.39 a is the Si-H bond length; 
and a ~ 3.820 A is the lattice constant. These structural parameters are taken from ab 
initio calculations [7]. The origin is placed at site A (figure 1). For the hydrogen atoms, 
we only consider up to first-nearest-neighbor interactions. We further assume that there 
is only a a bond between the Si and H atoms. Silicane has a hexagonal lattice with a 
basis composed of two Si atoms and two hydrogen atoms (figure 1). 

Our tight-binding model is constructed as follows: (i) the Si-Si TB parameters are 
taken from those of bulk Si {sp^ configuration); (ii) the Si-H TB parameters are obtained 
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Figure 1. Side view of the hydrogenated Si buckled sheet (silicane). Silicon atoms 
are located at sites A and B. Hydrogen atoms are located at sites C and D. The 
rectangle encloses the unit cell. 

Table 1. Tight-binding parameters for the hydrogenated Si sheet. Two-center 
Si-Si parameters are obtained from Grosso et al [20]. The H and Si-H parameters 
are obtained from Harrison [19] after applying the bond-length scaling for the Si-H 
parameter. All TB parameters are in eV. 



Si-Si 



S.(Si) 


-4.0497 


Ep 


1.0297 




-2.0662 




2.1184 


{ppa)tB 


3.1866 


{ppn)t^ 


-0.8867 




0.0000 


{spa)f^ 


0.0000 


{ppa)i^ 


0.8900 


ippTr)^^ 


-0.3612 


Si-H 


S.(H) 


13.6000 


(spa)^'^ 


6.21820 



by using the universal Harrison scaling rule [19] [iii) the two-center approximation 
is used to account for the slight distortion of the hydrogenated sheet from pure sp^ 
hybridization. A simple and fairly good parametrization of the bulk Si structure is that 
of Grosso et. al. [20]. These parameters reproduce fairly well the band structure of 
bulk Si, in particular in the KF direction, which is important for the Si nanosheets [3]. 
The full list of TB parameters is given in table 1. 

3. Results and Discussion 

We now present the results obtained from our tight-binding model. We first discuss the 
hydrogenated Si sheet. The hydrogenated Si sheet is an semiconductor with indirect 
gap Eg^ini{T — M) ~ 2.2 eV (figure 2). Our result is in between the 2.0 eV and 2.5 eV 
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Figure 2. Band structure of the hydrogenated Si sheet. The hydrogenated Si sheet 
is an semiconductor of indirect band gap -Eg.indir — M) ~ 2.2 eV. The Fermi level is 
set to zero (dashed line). 

reported previously [7, 17, 18]. The lowest point of the CB is at M and it is one- fold 
degenerate (figure 2 and table 2). The wave function at the M point of the lowest CB 
is given as follows (table 2): s (43%), (31%), p, (20%), and sh (6%). The highest 
point in the VB is at F and it is two- fold degenerate (figure 2 and table 2). The wave 
functions at the F point of the highest VB is pure px or pure Py (table 2). This is the 
same as was found for graphane [21]. 

The symmetries of the wave functions at F, M, and K for the hydrogenated Si 
sheet are given in table 2. Note that the K point {Esi.ii{K) = — 1.1 eV, 4.0 eV) is two- 
fold degenerate (table 2). This is expected since fully hydrogenating the Si sheet does 
not break the mirror symmetry of the honeycomb lattice. Another interesting result 
is that there is little hydrogen s orbital (sh) in the VB. Thus, the gap "opening " is, 
in reality, due to the hydrogen electrons filling up the originally empty lowest CB of 
silicene; i.e., the gap is mostly intrinsic to silicene. 

Next, we compare the band structures of graphane and silicane. The band 
structures of graphane and silicane differ significantly. First, graphane is an insulator of 
direct band gap [15, 16]. The gap occurs at F. Second, the band gap of graphane (5.4 eV) 
is far larger than that of silicane [16]. This is expected since the C-C bond length in 
graphane is greater than the Si-Si bond length of silicane [7, 15]. 

We now consider the band structure of Si-H NTs. Si-H NTs are fully characterized 
by the chiral vector (n, m), where n, m are integer numbers [22]. The band structure of 
Si-H NTs is obtained from that of silicane by quantizing the wave vector k along the 
chiral direction [22]. 

Armchair (n,n) and zig-zag (n,0) Si-H NTs with n even have concurrent and 
identical direct and indirect gaps equal to that of the hydrogenated Si sheet (2.2 eV) 
(figure 3(a) and (c)). We note that this dual behavior is different from the so-called 
mixed character previously reported [17] as the latter was for a nanosheet and was an 
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Table 2. Wave functions at F, M, and K for the hydrogenated Si sheet. The 
probabiUty densities of the electron states at the Si (H) atoms are given by the addition 
of the probability densities at sites A, B (C, D). 



i?si-H(r) [eV] Degeneracy Wave Functions 



16.3 


1 


z(14%),sj^(86%) 


16.0 


1 


z(ll%),s/f(89%) 


5.7 


2 


X (3%), y (97%); x {97%), y (3%) 


2.8 


1 


s(94%),z(5%),.s^f (1%) 





2 


X (100%); 2/ (100%) 


-1.5 


1 


s(4%),z(82%),sjy(14%) 


-5.0 


1 


s {6%), z (85%), SH {9%) 


-10.4 


1 


s(96%),z(4%) 


Es,.Ti{M) [eV] 


16.6 


1 


z {17%), SH (83%) 


16.5 


1 


z(15%),s/f(85%) 


7.9 


1 


y (100%) 


3.6 


1 


s (11%), X (70%), z (14%), (5%) 


2.2 


1 


s (43%), X (31%), z (20%), (6%) 


-0.8 


1 


s (12%), X (10%), z (66%), (12%) 


-1.4 


1 


2/(100%) 


-1.8 


1 


s (30%), X (1%), z (60%), SH (9%) 


-7.7 


1 


s(28%),x(67%),z(5%) 


-8.1 


1 


s(77%),x(19%),z(4%) 


Es,.n{K) [eV] 


16.6 


2 


z(17%),sff (83%) 






z(17%),Sff (83%) 


6.0 


1 


X (50%), 2/ (50%) 


4.0 


2 


s (19%), X (31%), y (31%), z (13%), sh (5%); 






s (19%), X (31%), y {31%,), z (13%), s// (5%) 


-1.1 


2 


s (16%), X (3%), y (3%), z (66%), (12%); 






s (16%), X (3%), y (3%), z (66%), s^/ (12%) 


-5.0 


1 


X (50%), 2/ (50%) 


-6.7 


2 


s (66%), X (15%), y (15%), z (4%); 



s (66%), X (15%), y (15%), z (4%) 



approximate phenomenon. Armchair (n,n) and zig-zag (t2,0) Si-H NTs with odd have 
direct gap equal to that of the hydrogenated Si sheet (figure 3 (b) and (d)). Therefore, 
Si-H NTs are semiconductors with either direct gaps, or with co-existing direct and 
indirect gaps of identical magnitude. 

The electronic behavior of Si-H NTs is explained as follows. We first consider 
armchair Si NTs. For armchair Si-H NTs {n,n), quantization of the wave vector k 
along the chiral direction gives their Brillouin zone (BZ) [22], 
27r o ^ 
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Figure 3. Band structure of (a) armchair (6,6), (b) armchair (7,7), (c) zig-zag 
(6,0), (d) zig-zag (7,0) hydrogenated Si nanotubes. Armchair (n,n) nanotubes with 
n even have simultaneously direct (F) and indirect (F X) gap equal to that of the 
hydrogenated Si sheet. Zig-zag (n,0) nanotubes with n even have simultaneously 
direct (F) and indirect (FX) gap equal to that of the hydrogenated Si sheet. For 
armchair nanotubes X= ±7r/a and for zig-zag nanotubes X= ±n/{\/3a). The Fermi 
level is set to zero (dashed line). The dotted line in (a) and (c) is a guide for the eye to 
show that the lowest conduction states at F and X are at the same energy. The dotted 
line in (b) and (d) is a guide for the eye to show that the lowest conduction states at 
F and X are not at the same energy. 



where, q = —n + 1, 0, and —ir/a < ky < n/a. Armchair {n,n) Si-H-NTs 
with n even have simultaneous and identical direct and indirect gaps since their BZ, 
equation (2), always crosses M and M' (figure 4 (a)), 

Armchair {n,n) Si-H NTs with n odd have only a direct gap since their BZ always 
crosses M but not M' (figure 4 (b)). 

We now consider zig-zag Si-H NTs. For zig-zag Si-H NTs {n,0), quantization of 
the wave vector k along the chiral direction gives their BZ [22], 

^ — kx H kjy, (3) 

where, q = — ra, 0, — 1 and — 7r/(-\/3a) < kx < Ti/{-\/^a). Zig-zag (n,0) Si-H NTs 
with n even have simultaneous and identical direct and indirect gaps since their BZ, 
equation (3), always crosses M, M' and M" (figure 4 (c)), 

, . 27r / 1 ^ 1 f \ , 27r / 1, ^ 1 \ . 27r 

Zig-zag (n,0) Si-H NTs with n odd have only direct gap since their BZ always crosses 
M" but not M, M' (figure 4 (d)). The above results on the gap dependence on chirality 
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(b) (d) 



Figure 4. Brillouin zone of (a) armchair (6,6), (b) armchair (7,7), (c) zig-zag (6,0), (d) 
zig-zag (7,0) hydrogenated Si nanotubes. Armchair (n,n) nanotubes with n even have 
simuhaneous direct and indirect gaps since its BriUouin zone crosses M and M'. Zig-zag 
(n,0) nanotubes with n even have simultaneous and identical direct and indirect gaps 
since its Brillouin zone crosses M, M' and M" . 

are much more complex than the Hamada rule [23] obeyed by non-hydrogenated carbon 
and silicon nanotubes [3]. 

We now compare our Si-H NTs results to those of Ref. [18]. Ref. [18] presents 
density functional tight-binding (DFTB) calculations of the Si-H NT gap. The DFTB 
calculation finds that Si-H NTs are semiconductors whose gap has little dependence on 
chirality. The gap of Si-H NTs rapidly converges to the gap of their hydrogenated Si 
sheet 2.5 eV). No reference is made as to whether the gap is direct or indirect. Even 
though the electronic behavior of our TB model is in fairly good agreement with that of 
DFTB, one important difference is that our TB calculation did not find a dependence 
of the size of the gap at all. We attribute this difference to curvature effects, which are 
neglected in our model. 
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4. Conclusions 

In summary, we performed a tight-binding calculation of the band structure of fully 
hydrogenated Si nanostructures. For the hydrogenated Si sheet (silicane), we find that 
it is a semiconductor with indirect band gap of about 2.2 eV. The indirect band gap 
occurs between the M and F points. Our study of the wave functions reveals that 
the gap is closely related to the band structure of silicene and is primarily due to a 
band-filling effect rather than to a gap opening. For Si-H NTs, we find that they are 
semiconductors with either a direct gap or with co-existing direct and indirect band gaps 
depending on chirality. Therefore, Si-H NTs do not follow Hamada's rule. This shows 
that hydrogenated Si nanotubes can have properties distinct from carbon nanotubes. 
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